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ABSTRACT 

Since  the  discovery  in  1924  that  we  live  in  a  separate  star  system,  now 
known  as  the  Milky  Way  Galaxy,  astronomers  have  worked  to  understand  the 
structure,  composition  and  dynamics  of  our  galaxy.  A  particular  puzzle  in 
recent  years  has  been  the  origin  and  dynamics  of  the  interstellar  gas  in  the 
halo  region  of  the  galaxy  high  above  the  visible  flat  stellar  disk.  The 
"chimney  model"  dominates  current  theories  of  the  source,  support  and 
ionization  of  the  halo  gas,  yet  no  one  has  actually  seen  a  galactic  chimney.  We 
selected  the  most  likely  candidate  for  a  galactic  chimney  in  our  galaxy  and 
observed  the  candidate  region  with  both  emission  and  absorption  spectra  over 
a  wide  wavelength  range.  We  were  granted  observation  time  to  observe  radio 
wavelength  emission  spectra  with  the  140  foot  National  Radio  Astronomy 
Observatory  telescope  at  Green  Bank,  West  Virginia,  and  ultraviolet  absorption 
spectra  with  the  National  Aeronautics  and  Space  Administration's 
International  Ultraviolet  Explorer  satellite.  These  original  observations, 
combined  with  previous  observations  in  radio,  optical  and  ultraviolet 
wavelengths,  provided  the  data  for  the  study.  The  investigation  takes  the 
novel  approach  of  observing  a  single  region  in  great  detail  at  many  different 
wavelengths  in  order  to  gain  understanding  of  the  composition,  ionization  and 
dynamics  of  the  halo  gas  in  the  region  of  a  possible  galactic  chimney. 
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I  Introduction 

Our  galaxy  is  a  spiral  galaxy  like  many  of  its  neighbors  (Figure  1).  The 
Sun  resides  about  8.5  kpc  (I  parsec  (pc)=  3.086x10^^  m  =  3.26  light  years)  from 
the  galactic  center  along  one  of  the  spiral  arms,  and  the  total  galactic  diame’er 
extends  approximately  50kpc.  The  galaxy  rotates  about  its  mass  centroid 
located  in.  the  dense  inner  galactic  nucleus.  A  nearly  spherical  mass  of  stars 


Figure  1:  Upper:  a  spiral  galaxy  seen  face-on  to  the  disk.  Lower:  a  spiral  galaxy  seen 
edge-on  to  the  disk. 
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called  the  bulge  surrounds  the  nucleus.  The  dense  region  in  the  plane  of  the 
galaxy  containing  the  spiral  arms  is  known  as  the  galactic  disk  while  the 
space  extending  far  above  and  below  the  disk  and  containing  only  a 
sprinkling  of  stars  is  called  the  galactic  halo. 

Vast  clouds  of  gas  and  dust  lie  between  the  stars  in  the  galactic  disk  and 
total  about  one  tenth  of  the  mass  of  the  galaxy.  Often  the  gas  is  diffuse  but 
occasionally  gas  is  also  found  in  dense  nebular  clouds,  such  as  the  Great 
Nebula  in  Orion,  that  shine  brilliantly  when  lit  by  nearby  stars.  One 
particularly  interesting  question  is  how  stars  and  interstellar  material  are 
related.  The  stars  in  the  Milky  Way  appear  to  have  grouped  by  age  with  the 
older  (Population  II)  stars  generally  found  in  the  halo  and  bulge  and  the 
younger  (Population  I)  stars  in  the  spiral  arms  of  the  disk.  Within  the  spiral 
arms,  the  interstellar  gas  and  dust  experience  the  pressure  required  to  form 
new  stars.  When  very  massive  stars  die  (run  out  of  effective  fuel  to  maintain 
fusion),  they  explode  violently  in  a  supernova.  With  tremendous  force,  the 
star  throws  out  the  gas  that  once  comprised  it  into  the  surrounding  area  of  the 
galaxy  where  the  gas  eventually  becomes  recycled  into  new  stars.  Because 
the  heavier  elements  only  form  within  stars  and  during  the  supernova 
explosion,  the  composition  and  metal  content  of  the  gas  tells  much  about  its 
history. 

While  most  of  the  interstellar  gas  appears  concentrated  in  the  disk, 
interstellar  clouds  of  both  neutral  and  ionized  gas  have  been  observed  high 
above  the  plane  of  the  galaxy  (Munch  and  Zirin  1961;  Albert  1983;  Albert  el.al. 
1992;  Danly  1989),  inspiring  many  questions  about  its  support  (what  keeps  the 
gas  up),  ionization  (the  degree  to  which  the  atoms  have  lost  electrons,  either 
from  absorbing  photons  or  from  collisions),  and  kinematics  (how  the  gas 
moves  around). 
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Two  major  theories  have  developed  explaining  the  support  and 
ionization  of  halo  gas.  In  the  "cold  halo"  theory,  cosmic  ray  pressure  provides 
the  support,  and  ionization  occurs  through  photoionization  (Chevalier  and 
Fransson  1984).  This  is  to  say  that  the  pressure  exerted  by  high  energy 
charged  particles  supports  the  gas  and  energy  imparted  by  the  light  from 
various  sources  ionizes  the  gas  while  no  energy  is  extracted  from  thermal 
means.  In  an  opposing  theory,  the  "hot  halo"  theory,  both  support  and 
ionization  are  achieved  by  thermal  means.  Heated  gas  rises  from  the  disk  into 
the  halo  and  then  cools  into  clouds  which  fall  back  into  the  disk  completing 
the  cycle  known  as  a  galactic  fountain  (Shapiro  and  Field  1975;  Brcgman  1980). 

The  galactic  chimney  model  explains  one  possible  outflow  mechanism  of 
the  fountain  model.  The  chimney  theory  proposes  that  correlated  groups  of 
.supemovae  in  the  disk  form  "superbubbles"  (Heiles  1986,  1987)  of  hot,  ionized 
gas,  which  may  expand  until  they  break  out  of  the  dense  disk  gas  into  the  halo, 
thus  forming  a  conduit  or  "chimney"  of  hot,  ionized  gas  between  the  disk  and 
the  halo.  The  gas  then  cools  and  rains  back  into  the  disk  as  in  other  fountain 
models  (Norman  and  Ikcuchi  1989).  Figure  2  illustrates  the  chimney  model. 


superbubble  chimney 


Figure  2;  A  qualitative  illustration  of  the  dynamic  cycle  of  a  galactic  chimney. 
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While  the  chimney  theory  has  been  well  developed,  no  one  has  actually 
seen  a  galactic  chimney.  We  explored  one  region  of  the  galaxy  where  a 
chimney  may  likely  be  located.  Two  methods  to  study  gas  in  the  halo  are  to 
observe  emission  spectra  of  neutral  hydrogen  and  to  observe  the  absorption 
spectra  of  other  ions  against  background  stars.  The  bulk  of  the  interstellar 
gas  is  neutral  hydrogen  (H  I),  which  is  observed  directly  in  emission,  giving  a 
"glow"  brighter  than  the  background  radiation  at  the  radio  wavelength  of  21 
centimeters.  Other  less  abundant  species,  such  as  calcium  and  titanium,  arc 
mixed  with  the  hydrogen  and  arc  observed  in  absorption.  Looking  at 
ab.sorption  spectra  is  somewhat  analogous  to  using  a  silhouette  (at  a  given 
wavelength)  against  the  light  from  a  distant  star  to  get  information  about  the 
matter  in  front  of  the  star.  Velocities  of  the  gas  are  determined  by  its  Doppler 
shift  where  the  ratio  of  the  change  in  the  wavelength  at  which  we  observe  a 
spectral  line  to  the  actual  wavelength  equals  the  ratio  of  the  velocity  of  the 
gas  to  the  speed  of  light  (Ak  /  X,  =  v/  3x10^  km/s).  Study  of  the  gas  at  many 
different  wavelengths  (and  therefore  many  different  energies)  reveals  much 
about  the  physics  of  the  gas,  giving  a  much  more  complete  picture  than  any 
one  wavelength  could  show. 

Cold  neutral  hydrogen  clouds  at  a  temperature  of  about  100  K  arc 
observed  with  radio  telescopes,  while  optical  telescopes  can  detect  trace 
element  ab.sorption  both  from  cold  and  from  warmer  ionized  gas  at  10,000  K. 
However,  because  of  their  atomic  structure,  most  chemical  elements  produce 
strong  absorption  in  the  ultraviolet  region  of  the  spectrum  and  can  be  studied 
over  a  wide  temperature  range  of  100  K  to  200,000  K  with  orbiting  ultraviolet 
satellites.  We  must  use  all  of  these  techniques  to  search  for  a  possible  galactic 
chimney. 
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II  The  Investigation 

Based  on  the  optical  and  ultraviolet  absorption  studies  of  Albert  and 
Blades  ct.  al.  (1992)  and  Danly  n989).  wc  selected  the  most  probable  location  of 
a  galactic  chimney  in  our  galaxy:  the  region  around  the  famous  halo  star 
HDl  19608.  A  number  of  unusual  characteristics  makes  this  the  best  region  for 
this  study.  This  is  the  only  area  to  show  a  marked  outflow  of  gas  in  low  ion 
(neutral  and  singly  ionized)  absoq)tion  spectra  and  in  H  1  emission.  At  the 
same  time  as  we  see  unusual  outflow  of  the  low  ions,  we  see  a  lack  of  the 
intermediate  velocity  (70  to  80  km/s)  infalling  low  ions  otherwise  found 
throughout  the  sky  (Danly  1989). 

In  addition  to  the  unusual  kinematics  of  the  gas  in  the  region,  the  line 
of  sight  toward  HDl  19608  also  has  unusual  ionization  characteristics.  HDl  19608 
is  the  only  high  latitude  star  to  show  absorption  from  the  very  highly  ionized 
ion  N  V,  an  indicator  of  hot  gas.  The  other  widely  studied  highly  ionized 
species,  C  IV  and  SI  IV,  arc  also  present. 

We  carefully  analyzed  absorption  data  toward  nine  stars  in  the  region, 
which  are  described  in  Table  1.  An  unusual  velocity  separation  between 
highly  ionized  species  and  low  ions  appears  in  the  absorption  spectra  toward 
HD119i.  ')8  and  two  other  stars,  HD121968  and  HD125924,  within  ten  degrees  of 
galactic  latitude  to  HDl  19608.  In  summary,  these  observations  suggest 
interesting  flows  of  both  hot  and  cold  gas,  which  may  be  indicative  of  a 
galactic  chimney.  Observations  in  other  regions  of  the  galaxy  reveal  none  of 


these  characteristics. 
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Table  1;  Star  Data 


Name 

MK  type 

1 

b 

ft  A 

Dec 

z 

V 

HD  1  17880 

B9 IV/V 

317 

43 

13  30  47.6 

- 18  15  24.9 

0.33 

9.1 

1 

HD  1  19608 

B1  lb 

320 

43 

13  41  48.2 

-[7  4\  10.8 

3.1 

7.56 

2 

HD  1  19786 

A2  V2 

322 

45 

13  42  53.2 

-15  31  01.3 

0.07 

6.2 

3 

HU123883 

B8/9  labp 

328 

41 

14  07  48.5 

-17  45  20.2 

8.7 

9.4 

2 

ljD*2  27  1  1 

02.5  V 

330 

62 

13  40  44.0 

♦0143  08 

2.5 

10.4 

3-7 

HD  1  18246 

l;3  me 

322 

- 1 

55 

13  33  66.7 

-05  54  03.6 

1 . 1 

8.07 

2 

HD  1  20086 

03  III 

330 

58 

13  4^  44  2 

-2  II  40  1 

0.85 

7.89 

8 

HD  123924 

02  IV 

338 

48 

14  20  37.7 

-8  01  15.9 

2.4 

9.68 

8 

HD  12 1968 

01 

334 

56 

4.3 

10.31 

8 

References:  (1)  Albert  1983.  (2)  Albert  et.al.  1992,  (3)  Tobin  1985,  (4)  Houk  1982 
(5)  Bidelman  1988.  (6)  Hill  1970.  (7)  Stetson  1985,  (8)  Danly  1989 


Table  1  lists  the  stars  included  in  the  study.  The  top  half  were  newly  observed  in  tlii.': 
Trident  project  while  the  lower  half  of  the  stars  were  observed  previously  by  Danly  (1989) 
and  Albert  et  al.(1992).  MK  type  is  the  spectral  classification  of  the  star  by  the  Morgan. 
Keenan  and  Kellman  classification  system.  I  and  b  are  galactic  longitude  and  galactic 
latitude  in  degrees.  Stellar  coordinates  arc  right  ascension  (RA)  in  hours,  minutes  and 
seconds  of  time  and  declination  (dec)  in  degrees,  minutes  and  seconds,  z  is  the  height  of 
the  star  above  the  galactic  plane  measured  in  kiloparsecs.  V  is  the  visual  magnitude  of  the 
s  I  a  r . 


The  region  of  study  extends  from  306°  to  332”  galactic  longitude  and 
+  34"  to  +64°  galactic  latitude.  Galactic  longitude  (I)  is  measured  starting  with 
the  position  of  the  sun  relative  to  the  center  of  the  galaxy  at  0°  and 
continuing  around  the  disk  counterclockwise  from  north  to  360°.  Galactic 
latitude  (b)  is  measured  going  from  0°  to  90°  above  the  plane  of  the  galaxy  and 
()"  to  -90“  below  the  galactic  plane.  Figure  3  shows  the  halo  stars  in  the  region 
as  they  appear  relative  to  one  another  on  the  sky  measured  in  the  celestial 
coordinates  of  right  ascension  and  declination. 
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Figure  3:  Background  stars  observed  in  the  region  plotted  in  coordinates  of  right 
ascension  (hours  and  minutes)  and  declination  (degrees)  as  they  are  seen  in  the  sky  in  the 
constellation  Virgo. 


Figure  4  portrays  the  galaxy  from  a  side  view  looking  edge  on  to  the 
disk  as  in  the  bottom  part  of  the  galaxy  diagram  in  Figure  1,  which  shows  the 
relative  distances  of  the  stars  above  the  disk.  In  order  to  obtain  interstellar 
absorption  lines  we  must  use  a  suitable  star  as  a  background  source.  Hot  stars 
have  an  intense  continuum  and  few  absorption  lines  in  the  atmosphere  of  the 
star  it.sclf  which  might  cause  confusion  with  the  interstellar  lines.  Also 
important  in  the  selection  of  background  stars  is  choosing  distant  stars  of 
sufficient  brightness  to  make  observations.  In  order  to  establish  constraints 
on  the  distance  to  the  gas,  we  observed  nearby  foreground  stars  (HDl  19786 
and  HD  117880)  and  a  star  more  distant  than  HDl  19608  (HD123884). 


Figure  4;  The  stars  in  the  region  toward  HDl  19608  in  the  halo  of  the  galaxy  with  vertical 

height,  z,  in  kiloparsecs  above  the  plane.  Hie  angle  of  the  stars  above  the  plane  is  the 

galactic  latitude.  (Galactic  disk  is  not  to  scale). 

If  an  absorption  feature  is  observed  superimposed  in  the  spectrum  of  a 
distant  star  but  not  toward  a  nearby  star  in  the  same  direction,  then  wc  know 
that  the  gas  of  the  feature  must  be  located  beyond  the  foreground  star  and  in 

front  of  the  distant  star,  unless  the  gas  is  extremely  patchy  and  varies  in  the 

very  small  angle  between  the  stars  in  the  line  of  sight.  In  order  to  establish  a 
vertical  limit  of  the  feature  (how  high  the  feature  extends  in  galactic  latitude) 
wc  observed  the  high  latitude  star  BD+2  2711  at  62*  ‘ 

The  purpose  of  this  investigation  is  to  probe  the  region  suspected  to  be  a 
possible  galactic  chimney  to  determine  spatial  constraints  (localize  the 
unusual  feature  in  space),  obtain  ionization  and  abundance  information,  and 
observe  kincmajc  effects.  The  next  three  sections  will  discuss  and  analyze 
observations  of  the  region  in  radio,  optical  and  ultraviolet  wavelengths.  We 
then  compare  the  results  to  the  chimney  model  to  support  or  refute  the  theory 


Ill  Radio  Observations  and  Analysis 


Observations  of  neutral  hydrogen  emission  give  a  spatial  picture  of  the 
location  of  interstellar  gas  on  the  plane  of  the  sky  because  emission  can  be 
observed  in  any  direction,  not  just  in  lines  of  sight  where  background  stars 
are  available.  The  energy  transition  which  results  in  the  emission  of  a  photon 
with  a  radio  wavelength  of  21  centimeters  from  a  hydrogen  atom  occurs  when 
the  electron  in  the  first  orbital  (n=l)  and  with  angular  momentum  of  0  (1=0) 
changes  its  spin  orientation  with  respect  to  the  spin  of  the  proton  in  the 
nucleus  (Figure  5). 


F^l 


Figure  5:  The  energy  level  diagram  of  the  ground  state  of  neutral  hydrogen. 


The  electron  with  its  spin  parallel  to  the  proton  (total  spin  quantum  number 
(F)  =  electron  spin  +  proton  spin  =1/2  +  1/2  =  1)  "flips"  to  spin  antiparallel  (F  = 
1/2  -  1/2  =  0).  The  energy  change  is  very  small,  only  5.9x10'^  eV,  resulting  in 
a  very  long  wavelength  of  emission.  21  cm.  This  example  of  hyperfine 
structure  (evidence  of  electron-nucleus  spin  interaction)  occurs  very  rarely, 
with  a  probability  of  transition  of  one  time  in  10^  years.  Therefore,  the  H  1 
hyperfine  transition  is  considered  "forbidden"  in  the  relatively  dense 
laboratory  conditions  on  earth.  In  a  dense  gas,  the  spin  flip  is  caused  by 
collisions  long  before  the  low-probability  spontaneous  transition  can  occur. 


In  interstellar  space,  however,  the  average  density  of  one  atom/cm^  is  so  low 
that  the  atoms  remain  in  the  excited  state  and  the  "forbidden"  transition  is 
common.  Hydrogen  makes  up  over  70  percent  of  all  matter,  and  the  enormous 
volume  containing  this  most  abundant  clement  makes  the  spectral  line  very 
strong  and  easy  to  measure. 

We  constructed  contour  maps  of  the  emission  of  neutral  hydrogen  over 
the  northern  galactic  sky  from  the  all  sky  survey  done  by  Bell  Labs  (Stark 
1992).  The  only  marked  region  of  outflow  above  the  galactic  plane  was 
strongest  at  a  velocity  of  25  km/s  and  located  in  the  region  of  our  background 
stars.  Figure  6  shows  the  contour  map  of  the  emission  strength  of  the  25km/s 
outflowing  gas  over  a  polar  view  of  the  northern  galactic  sky. 
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Figure  6:  Bell  Labs  survey  of  the  2lcm  radio  emission  in  the  northern  galactic  sky  mapped 
in  contours  indicating  the  25  km/s  outflow  strength.  The  arrow  indicates  the  outflow 
region  where  MDl  19608  is  located. 


Observations 
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The  Bell  Labs  survey  has  the  low  spatial  resolution  of  three  degrees,  so 
in  order  to  fully  explore  the  spatial  characteristic  of  the  outflow  region  at 
high  resolution,  we  applied  for  observing  time  for  this  Trident  project  on  the 
140  foot  radio  telescope  at  ihe  National  Radio  Astronomy  Observatory  (NRAO)  in 
Green  Bank,  West  Viiginia,  which  is  one  of  two  radio  observatories  supported 
by  the  National  Science  Foundation.  Completed  in  1965,  the  instrument  used 
for  observation  is  the  largest  equatorial  mount  telescope  in  the  world  (Figure 
7).  NRAO  granted  us  time  on  three  occasions  and  we  conducted  observations  on 
1-2  April  1991,  10-12  August  1991,  and  20-21  January  1992. 


Figure  7:  The  140  fool  telescope  of  the  National  Radio  Astronomy  Observatory  in  Green 
Bank,  West  Virginia. 


The  telescope  focuses  the  radio  waves  using  a  solid-surface,  140  foot 
diameter  paraboloid  at  a  focal  ratio  of  0.43  onto  a  receiver  about  the  size  of  a 
refrigerator.  Two  antennas  receive  the  signal  at  perpendicular  polarizations. 
A  complex  system  of  cryogenics  maintains  the  internal  electrical  system  noise 
below  the  cosmic  background  radiation  of  3  degrees  Kelvin.  The  150  ton 
telescope  dish  rotates  on  a  single  very  large  ball  bearing  which  is  kept 
constantly  lubricated  with  a  thin  layer  of  oil  on  which  the  dish  "floats".  The 
telescope  has  a  pointing  accuracy  and  repeatability  of  10  arcscconds  (  .003 
degrees).  Positioning  parameters  and  information  are  written  into  a  computer 
format  which  mimics  the  original  punched  cards  first  used  at  the  telescope. 

The  telescope  operates  automatically  according  to  computer  input,  and  a 
professional  telescope  operator  monitors  the  computer  and  engineering 
systems  and  may  enter  commands  manually.  The  observer  directs  the 
observations  and  can  view  raw  spectra  on  a  computer  work  station  within 
minutes  of  the  time  they  are  taken. 

We  utilized  our  observing  time  to  scan  the  region  around  HDl  19608, 
scanning  along  the  horizontal  direction  across  the  sky  in  right  ascension. 
Spectra  were  taken  at  a  rate  of  12  per  minute  with  a  20  areminute  beam  width. 
The  right  ascension  "strips"  with  one  spectrum  every  10  areminutes  were 
taken  at  30  areminute  intervals  in  the  vertical  direction  of  declination  to 
construct  a  map  of  the  radio  spectra  at  the  high  spatial  resolution  of  20 
areminutes.  The  initial  survey  indicated  a  strong  feature  of  outflowing  gas  at 
25  km/s,  and  the  subsequent  observations  extended  the  map  to  determine 
boundaries  of  the  flow  feature  and  to  explore  the  higher  latitudes  around  the 


star  BD+2  2711. 


As  with  all  radio  telescopes,  the  140  foot  has  sidelobe  radiation  beams  in 
addition  to  the  main  beam  pattern.  In  halo  gas  in  which  the  hydrogen 
emission  is  weak  compared  to  that  of  the  disk,  accidental  measurement  of 
scattered  signal  picked  up  in  the  sidelobes  can  seriously  degrade  the  data. 
Therefore,  the  program  included  observations  toward  each  star  to  measure 
the  stray  radiation  emitted  into  the  telescope  beam  sidelobes  and  to  compare 
these  measurements  with  standard  spectra.  Comparisons  to  standard 
measurements  indicate  no  significant  stray  radiation  detections  in  our 
observations. 


Radio  Analysis 

We  observed  over  3000  individual  positions  (times  2  polarizations  gives 
over  6000  spectra)  toward  the  region.  The  first  step  in  analysis  of  the  H  1 
spectra  is  to  determine  a  baseline  which  sets  the  background  continuum  to  a 
value  of  zero.  Almost  the  entire  region  has  a  very  flat  background  radio 
continuum,  so  we  used  primarily  a  DC  offset  to  simply  move  the  base  of  the 

spectra  to  a  zero  flux  value  (Figure  8).  The  Unipops  radio  data  analysis  software 

package  (written  by  Ronald  Maddalena  of  NRAO)  used  to  baseline  the  data 
allows  baseline  fits  to  polynomials  in  cases  where  the  background  is  not  flat. 

Parameters  specified  in  the  baselining  procedure  are  the  beginning  and  end 
velocities  to  include  in  the  fit  and  the  width  of  the  regions  to  be  averaged  to 
dcicrminc  the  baseline  offset.  Figure  8  shows  an  example  spectrum  before  and 

after  baselining  with  a  DC  offset.  The  boxed  regions  indicate  the  points  over 
which  the  average  value  for  the  offset  was  computed.  The  remainder  of  the 
data  analysis  was  accompli.shcd  on  the  Hl-IDL  software  package  (written  by  K. 
D.  Kimtz  of  the  Space  Telescope  Science  Institute).  The  ba.sclincd  spectra  arc 


transferred  from  the  magnetic  tape  on  which  they  are  stored  at  the  telescope 
in  ASCII  format  into  a  disk  file  of  a  format  readable  by  the  HI-IDL  software 
package. 
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Figure  8:  An  example  .spectrum  before  baselining  (at  lop)  and  after  baselining  the 
continuum  (on  bottom).  The  horizontal  scale  is  velocity  in  km/s  and  the  vertical  scale  is 
the  actual  number  of  neutral  hydrogen  aioms/cm^  (column  density). 


To  convert  from  measured  antenna  strength  to  brightness  temperature 
(the  standard  unit  of  radio  emission)  and  to  correct  for  absorption  in  the 
atmosphere,  we  use  the  following  relation  (Williams  1972): 

Tb  =  f  Ta 

T  =  0.0101/cos  (Zenith  angle) 
f  =  1.0818  for  antenna  1  (scan  numbers  ending  in  .01) 
f  =  1.0515  for  antenna  2  (scan  numbers  ending  in  .02) 

T  b  is  brightness  temperature  on  the  Berkeley  scale 
T  a  is  antenna  strength  (vertical  axis  of  raw  spectral  plots) 

The  value  of  f  is  a  function  of  the  beam  pattern  of  the  particular  telescope  and 
antenna  used  and  has  been  determined  experimentally  for  the  Green  Bank  140 
foot  telescope.  The  atmospheric  absorption  is  a  small  correction  at  radio 
wavelengths  and  depends  on  zenith  angle  (the  angle  from  overhead  to  the 
horizon).  The  correction  will  be  greater  in  observations  at  a  large  zenith 
angle  than  those  high  overhead  since  we  are  looking  through  more  of  the 
atmosphere  at  the  horizon  (Figure  9). 


Figure  9;  Qualitative  demonstration  that  the  path  length  of  a  wave  through  the  atmosphere 
(indieaicd  by  the  dotted  region)  is  greater  near  the  horizon. 


Wc  then  divided  the  region  into  small  areas  or  bins.  The  ideal  lower 
limit  of  the  bin  size  is  the  beam  width  of  the  telescope  but  the  practical  lower 
limit  was  the  need  to  include  more  than  one  spectrum  per  bin.  Fortunately, 
our  beam  size  allowed  an  average  of  three  spectra  per  bin  (Figure  10).  We 
divided  the  velocities  into  segments  to  look  at  the  strength  of  the  emission  in  a 
discrete  velocity  increment.  Shown  below  on  the  left  are  the  position  of  the 
spectra  indicated  by  dots  and  on  the  right  the  bin  grid  overlap  indicating  the 
spectra  in  each  spatial  bin. 
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Figure  10:  Left:  The  position  of  individual  scans  as  dots.  Right:  The  grid  overlay  showing 
the  spectra  in  each  spatial  bin. 


Using  the  two  dimensional  array  called  a  channel  map,  which  contains 
the  values  of  the  integral  of  the  brightness  temperature  over  the  given 
velocity  increment  Av  and  bin  space  x  and  y,  we  get  column  density  (N)  at 
each  spatial  bin; 


v^Av 

N  =  1.823  X  X  X  X  y  X  j  T(v)dv 

The  column  density  is  the  actual  number  of  atoms  along  the  line  of  sight 
producing  the  observed  spectral  feature.  Contour  maps  of  the  column  density 
in  the  observed  region  were  constructed  for  each  5  km/s  increment  from  -100 
km/s  to  100  km/s. 

Of  particular  significance  is  the  strong  outflow  apparent  at  the  25.5 
km/s  (A  velocity  of  2.5  km/s  on  either  side)  outflow  map  (  Figure  11).  Black 
indicates  no  observation  was  taken  in  the  region  or  no  gas  was  seen  in  that 
velocity  range.  The  white  square  outlines  the  single  bin  space  in  the  observed 
region  containing  no  spectra.  White  is  the  lowest  level  of  column  density, 
wliich  is  at  the  lower  limit  of  the  signal  to  noise  detection  limit.  The  levels  are 
filled  in  successively  darker  shades  with  each  column  density  level  as  shown 
in  the  legend. 

The  position  of  several  of  the  observed  stars  are  shown  on  the  contour 
diagram.  The  map  shows  the  strongest  outflow  (darkest  level)  in  the  area  of 
HDl  19608  with  weak  oulflow  on  the  lower  longitude  side  at  the  right  around 
the  position  of  the  star  HD  123884.  At  this  velocity,  the  column  density 
generally  decreases  as  the  latitude  increases  to  the  star  BD+2  2711.  Included  in 
.Appcndi,\  A  arc  all  of  the  velocity  contour  maps  over  the  velocities  from 
-100  km/s  to  -t-100  km/s.  The  contour  levels  arc  the  same  for  all  maps. 
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STAR  LEGEND 

■  HD  119608  •  HD  1 23084  □  HD  I  19786  O  HD  11  7880  /^IbD-2  2711 


Figure  1 1 :  The  contour  map  of  neutral  hydrogen  emission  strength  at  the  strongest 
positive  velocity  feature  at  the  veloeity  of  25  km/s. 


20 


Figure  12  shows  the  clean  "on  star"  spectra  toward  each  of  the  observed 
stars  (Albert  ct  al.  1992).  The  25  km/s  feature  can  be  easily  distinguished  in 
the  emission  spectra  toward  the  stars  HDl  19608,  HDl  17880  and  HDl  19786.  The 
25  km/s  emission  is  very  weak  toward  HDl 23884  and  relatively  weak  toward 


Figure  12:  The  individual  neutral  hydrogen  emission  spectra  from  the  direction  of 
HD119608,  HD117880,  HD123884,  BD+2  2711.  and  HD119786.  The  spectra  each  cover 
.slightly  different  ranges  of  velocity  as  shown  on  the  horizontal  scale  in  km/s.  The 
vertical  scale  is  brightness  temperature. 


IV  OPTICAL  OBSERVATIONS  AND  ANALYSIS 
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While  the  21  cm  neutral  hydrogen  emission  gives  a  two  dimensional 
map  of  the  velocities  of  hydrogen  gas  ir  space,  it  does  not  tell  us  how  far 
above  the  disk  the  interstellar  gas  clouds  are  located.  In  order  to  determine 
distance  limits,  we  use  very  high  velocity  resolution  optical  absorption 
spectra.  The  K  line  of  singly  ionized  calcium  at  3933.633  A  toward  HDl  19608, 

HDl  17880,  HDl  18246,  HD  119786.  and  HD123886  and  the  line  of  singly  ionized 
titanium  ai  3383.761  A  toward  the  first  four  of  these  stars  were  observed  in  a 
large  survey  of  halo  gas  by  Albert  and  Blades  (1992).  The  optical  absorption 
observations  were  conducted  using  a  charge  coupled  device  (CCD)  detector 
with  the  cchelle  grating  of  the  coude  spectrograph  of  the  Canada  France 
Hawaii  Telescope  on  Mauna  Kca,  Hawaii.  The  resulting  spectra  achieved  an 
average  velocity  resolution  of  5  km/s.  Figure  13  shows  the  singly  ionized 
calcium  (Ca  II)  line  toward  each  of  the  stars  and  the  line  of  singly  ionized 
titanium  (Ti  II)  toward  HDl  19608. 

The  digitized  spectra  arc  shown  as  plots  of  intensity  versus  velocity.  In 
most  interstellar  studies,  velocity  is  measured  with  respect  to  the  Local 
.Standard  of  Rest  (LSR).  a  reference  frame  defined  by  averaging  the  motions  of 
nearby  stars.  As  gas  moves  away  from  us  out  of  the  plane  of  the  galaxy,  its 
spectral  lines  arc  red  shifted  to  longer  wavelengths  and  positive  LSR 
velocities.  Infalling  gas  will  show  blue  shifted  absorption  with  shorter 
wavelengths  and  negative  LSR  velocities.  The  horizontal  axis  of  the  spectra  is 
the  velocity  of  this  Doppler  shift  in  km/s.  The  vertical  axis  of  an  absorption 
line  shows  the  intensity  of  the  detected  light,  with  the  background  continuum 
starlight  normalized  to  one.  The  line  "drops"  in  intensity  at  velocities  where 
the  gas  has  absorbed  some  of  the  background  light. 
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Figure  13:  The  optical  absorption  lines  of  singly  ionized  calcium  towards  stars  at  varying 
heights  above  the  disk  and  the  spectrum  of  singly  ionized  titanium  toward  HDl  19608. 


Singly  ioni/xd  calcium  i.s  one  of  the  strongest  and  best  studied  of  the 
optical  interstellar  species.  Ti  11  is  mor",  difficult  to  observe,  as  detector 
sensitivities  fall  off  considerably  at  this  wavelength,  so  titanium  can  only  be 


observed  in  the  brightest  stars. 
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Optical  Analysis 

Analysis  of  the  absorption  spectra  gives  us  information  about  distance 
and  ionization  of  the  clouds.  It  is  clear  from  the  spectra  shown  in  Figure  ’ 
that  very  little  gaseous  calcium  exists  up  to  a  distance  of  70  pc  above  >c 
galactic  plane,  and  a  weak  zero  velocity  feature  is  located  between  70  and 
330  pc.  A  very  interesting  outflowing  cloud  centered  at  +25  km/s  appears 
between  330  pc  and  3.1  kpc,  which  is  the  distance  to  HDl  19608.  It  is  also 
ob.scrvcd  as  a  weak  feature  in  the  spectrum  of  Ca  II  toward  a  neighboring  star, 
HDl  18246,  at  a  height  of  1.1  kpc  above  the  galactic  plane  (Albert  1992).  The  25 
km/s  feature  beeomes  stronger  in  the  spectrum  of  the  more  distant  star 
HD  123884,  showing  that  more  outflowing  gas  lies  between  HDl  19608  and 
HD123884. 

The  observed  features  in  the  spectra  are  actually  the  superposition  of 
absorption  from  many  clouds  at  varying  distances  and  perhaps  with  varying 
physical  conditions  along  the  line  of  sight.  To  understand  the  properties  of 
the  different  clouds  that  form  an  absorption  line  requires  the  detailed 
analysis  of  the  line  profile  to  separate  individual  clouds.  Cloud  modelling  can 
only  be  done  with  high  resolution  data. 

The  modelling  was  accomplished  using  the  profile  fitting  routine  of  the 
KDAF  spectral  data  reduction  software  package  written  by  K.D.  Kuntz  of  the 
Space  Telescope  Science  Institute.  The  routine  takes  as  input  a  model  of  many 
clouds  along  the  line  of  sight  each  specified  by  central  cloud  velocity,  Doppler 
broadening  (which  depends  on  cloud  temperature),  and  column  density.  The 

routine  uses  the  cloud  parameters  to  create  the  line  profile  of  a  given  spectral 
line  as  it  niighi  be  observed  with  a  particular  instrument.  This  simulated  line 


profile  can  be  compared  to  the  measured  line  profile.  To  construct  a  model,  we 
varied  the  parameters  to  get  the  best  possible  fit  with  only  one  cloud,  then  two, 
etc.,  adding  clouds  until  the  modelled  profile  matched  the  observed  profile  to 
within  the  signal  to  noise  limit  of  the  data.  The  resulting  models  of  HDl  19608 
and  HD123884  each  contained  eight  clouds.  The  positions  and  Doppler  widths  of 
the  clouds  remained  the  same  for  both  stars,  but  the  column  densities  were 
increased  for  HD  123884  to  account  for  the  greater  distance  of  that  star.  The 
profile  fitting  technique  gave  us  column  densities  of  the  calcium  at  the 
distances  to  each  star  and  showed  outflowing  calcium  clouds  at  velocities 
around  25  km/s,  high  above  the  plane  of  the  galaxy. 

To  understand  the  physical  conditions  in  each  cloud,  it  is  necessary  to 
determine  whether  the  hydrogen  gas,  which  forms  the  bulk  of  the  interstellar 
material,  is  neutral  or  ionized.  Observations  of  Ca  II  determine  neither  the 
gaseous  abundance  of  calcium  nor  the  level  of  ionization  in  the  cloud,  since 
both  Ca  II  and  the  unobservable  Ca  III  may  coexist  in  neutral  or  ionized 
hydrogen  regions.  Singly  ionized  titanium,  however,  is  a  unique  tracer  of 
neutral  hydrogen  due  to  the  almost  exact  coincidence  of  the  ionization 
potentials  of  Ti  II  and  H  I.  A  cloud  observed  in  Ti  II  is  almost  certainly  a 
neutral  hydrogen  region  and  the  strength  of  the  Ti  II  absorption  line  docs 
reflect  its  true  gaseous  abundance.  The  Ti  II  line  toward  HDl  19608  matches 
wcl'  with  the  Ca  II  line  and  also  shows  the  25  km/s  outflow  feature  indicating 
that  both  species  are  indeed  in  the  same  cloud  as  the  neutral  hydrogen.  On 
the  other  hand,  the  neutral  hydrogen  spectrum  observed  in  radio  emission 
toward  HD123884  (which  shows  the  strongest  outflow  of  Ca  11)  is  surprisingly 
weak  in  the  positive  velocities.  Although  HD  123884  is  too  faint  to  obtain  Ti  II 
spectra  with  current  technology,  the  lack  of  H  I  indicates  that  the  line  of 
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sight  toward  HDl 23884  does  not  go  through  the  neutral  cloud  itself  but  through 
the  ionized  outer  edge.  This  is  a  very  interesting  and  unique  opportunity  to 
probe  the  detailed  inner  and  outer  structure  of  an  interstellar  cloud. 

Titanium  and  calcium  also  share  the  useful  property  of  being  refractory 
elements.  Such  elements  can  very  easily  condense  out  of  the  interstellar  gas 
onto  solid  dust  grains  and  are  easily  released  back  into  the  gas  phase  if  the 
grains  are  destroyed  by  a  passing  shock  front  (Jenkins  1986).  Absorption 
lines  of  Ca  II  and  Ti  II  would  therefore  be  enhanced  in  gas  that  had  been 
shocked  by  a  supernova  remnant  or  had  been  propelled  through  high  velocity 
flows  such  as  a  galactic  fountain  or  chimney.  The  lines  observed  toward  the 
stars  HDl  19608  and  HD123884  are  strong  compared  to  typical  disk  gas  and  arc 
consistent  with  a  chimney  scenario. 

V  Ultraviolet  Observations  and  Analysis 

The  optical  observations  were  conducted  with  the  most  advanced  ground 
based  technology  at  high  resolution  which  allows  detailed  analysis  of 
individual  velocity  components  in  the  absorption  spectra.  The  high 
resolution  is  critical  to  the  identification  of  the  weak,  narrow  components  of 
numerous  clouds  which  arc  superimposed  in  a  typical  interstellar  absorption 
line  and  to  distinguish  the  lines  of  many  different  elements  found  in  a  very 
small  range  of  wavelength.  However,  only  a  very  limited  number  of  ions 

ah.sorb  in  the  wavelengths  visible  to  the  human  eye.  Due  to  atomic  structure, 
many  more  ions  have  strong  resonance  ab.sorption  lines  in  the  ultraviolet 
hand  of  the  spectrum.  Therefore,  even  though  ultraviolet  spectra  (with 
average  resolution  of  25  km/s)  cannot  achieve  the  high  resolution  of  optical 
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observations  (5  km/s),  the  spectra  in  the  ultraviolet  wavelengths  from  1190  to 
2070  A  proved  essential  in  providing  ionization  and  abundance  information. 

As  part  of  this  prohect,  we  applied  to  NASA  in  January  of  1991  for  use  of 
satellite  time,  and  were  able  to  obtain  spectra  with  the  International 
Ultraviolet  Explorer  (lUE)  satellite  in  three  shifts  from  26  to  28  June  1991.  We 
observed  HD123884  in  two  shifts  and  BD+2  2711  in  the  third  shift  from  the 
ground  station  at  the  NASA  Goddard  Space  Flight  Center  in  Greenbcll, 
Maryland. 

Launched  in  1978.  the  lUE  satellite  orbits  the  earth  in  an  elliptical 
geosynchronous  orbit  in  which  the  craft  is  always  visible  from  one  of  two 
ground  stations,  the  Goddard  Space  Flight  Center  operated  by  NASA  or  the 
Villafranca  Satellite  Tracking  Station  near  Madrid,  Spain  operated  by  the 
European  Space  Agency.  Goddard  controls  the  satellite  during  two  eight  hour 
shifts  each  day  and  VILSPA  operates  the  third  shift. 

Two  large  solar  panels  power  the  spacecraft  with  a  backup  battery  for 
periods  when  the  satellite  travels  through  the  earth’s  shadow  (Figure  14).  The 
on-board  computer  controls  telescope  pointing,  execution  of  movements,  and 
cxpo.surc  time.  The  satellite  has  three  axis  stability  (roll,  pitch,  and  yaw)  and  a 
nominal  pointing  accuracy  of  one  arcsecond.  Attitude  control  is  accomplished 
using  the  two  remaining  operable  gyroscopes  of  the  original  six,  together 
with  the  Fine  Sun  Sensor  which  detects  the  position  of  the  sun  (Sonneborn  et 
al.  1987).  Reaction  wheels  are  used  to  store  angular  momentum  and 
accomplish  spacecraft  rotation.  Therefore,  after  a  movement  of  the  telescope 
into  position  (called  a  slew),  and  before  observations  are  made,  the  wheels 
often  must  be  "unloaded",  using  small  hydrozene  jets,  to  establish  proper 
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Figure  14:  The  International  Ultraviolet  Explorer  Satellite  operated  by  NASA  and  ESA. 


During  observations,  a  telescope  operator  and  resident  astronomer 
control  the  satellite  and  provide  technical  assistance,  with  directions  from  the 
guest  observers.  Operations  Control  Center  personnel  monitor  all  of  the 
engineering  functions  of  the  satellite  in  a  building  nearby  the  Telescope 
Operations  Center.  The  ground  station  transmits  commands  to  the  satellite  at 
VHF  (Very  High  Frequencies),  while  the  satellite  in  turn  relays  information  at 
a  rate  of  20  kilobytes/scc  in  the  S-band  frequencies.  An  18  meter  antenna  at 
the  International  Ultraviolet  Explorer  (lUE)  tracking  station  on  Wallops 
Island.  Virginia,  receives  the  information  from  the  satellite  and  forwards  the 
data  to  Goddard  via  a  communications  satellite. 

The  lUE  spacecraft  has  four  cameras  and  two  spectrographs.  The  two 
spectrographs,  one  for  short  and  the  other  for  long  ultraviolet  wavelengths, 
and  each  has  a  primary  and  a  redundant  camera.  We  m.adc  each  of  our 
observations  using  ihc  SWP  (Short  Wavelength  Prime)  camera.  The  camera 
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electronically  counts  the  number  of  photons  in  each  wavelength  registered 
in  photomultiplier  tubes.  Over  the  exposure,  the  image  is  integrated  in  the  SEC 
Vidicon  section  of  the  camera.  The  camera  retains  its  image  until  the 
information  is  read  down  to  the  tracking  station  using  a  raster  scan  of  768x768 
pixels.  The  video  signal  is  digitized  into  256  discrete  levels  (0  to  255  Data 
Numbers  (DN))  by  an  eight  bit  analogue  to  digital  converter.  Real  time 
transmission  occurs  to  the  ground  station  and  a  residual  image  remains  on 
the  camera  until  after  a  PREP  is  executed.  A  PREP  consists  of  exposing  the 
image  with  a  tungsten  flood  lamp  to  get  a  reproducible  electronic  pedestal  of 
15  to  40  DN.  To  check  the  camera,  a  reading  is  then  taken  with  the  defocused 
beam  and  the  exposure  and  reading  are  repeated.  The  PREP  procedure 
basically  exposes  the  image  to  a  constant  starting  background  to  get  rid  of  the 
leftover  image  from  the  previous  observation.  During  observations, 
incoming  light  is  focused  onto  an  echelle  spectrograph  which  separates  the 
light  into  its  various  wavelengths  giving  a  series  of  strips  called  orders. 

Photons  at  each  wavelength  then  register  in  photomultiplier  tubes.  Careful 
monitoring  of  background  radiation  is  required  throughout  an  observation  to 
avoid  overexposing  and  damaging  the  instrument. 

Ultraviolet  Observations 

To  conduct  an  observation  toward  a  star,  we  make  a  finding  chart  from 
the  Palomar  Sky  Survey  prints  to  identify  the  star  in  the  image  from  the 
satellite's  Fine  Error  Sensor  (FES).  The  FES  uses  a  photocathode  that  has  a  peak 
response  in  the  visible  wavelengths  to  provide  a  television  image  of  the  field 
of  view.  When  the  target  has  been  properly  identified,  the  exposure  is  begun. 
At  the  end  of  the  exposure  the  image  is  read  down  and  can  be  viewed  in  the 
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form  of  a  phoiowriie  (Figure  15).  The  digital  spectrum  is  stored  from  the 
computer  onto  magnetic  tape  for  future  analysis.  Two  observations  were  taken 
of  HD123884;  a  420  minute  exposure  on  26  June  1991  recorded  in  the  file 
"SWP41921hlg"  and  a  430  minute  exposure  on  28  June  1991  in  "SWP41944hlg". 
One  320  minute  exposure  of  BD+2  2711  is  recorded  in  "SWP41932hlg".  All  of  the 
observations  made  on  lUE  become  public  record  in  .archives  a  year  after  the 
ob.scrvation.  HDl  19608  wa.s  observed  in  1979  and  we  used  ihe  archived  spectra 
toward  this  star  in  addition  to  our  own  observations  of  HD123884  and  BD+2  2711 
for  our  ultraviolet  analysis. 


lUESOC  9843  FILE;1  SWP41921  HD123884  GPiNLD 
RC:31  6/27/91  8:17  fiM 


Figure  15;  A  phoiowritc,  which  is  a  photographic  representation  of  the  ob.served  spectrum 
of  the  .star  HD  123884.  Each  strip  is  an  order.  Evenly  spaced  white  marks  are  reseaux 
(alignment  marks).  Black  spots  arc  cosmic  ray  hits. 
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Ultraviolet  Analysis 

The  spectral  files  from  the  magnetic  tape  of  the  observing  run  were 
transferred  onto  disk  files  on  the  "Sol”  computer  system  at  The  Space  Telescope 
Science  Institute  where  all  of  the  analysis  of  the  ultraviolet  spectra  was 

accomplished  using  the  KDAF  spectral  data  reduction  software  package.  The 
analysis  of  the  ultraviolet  spectra  requires  a  series  of  corrections  and 
procedures  in  order  to  obtain  each  individual  interstellar  absorption  line  for 
study.  We  extracted  each  spectral  line  of  interest  as  specified  in  an  atomic  line 

data  file  [Sec  Appendix  B:  Atomic  Line  Data  File]  and  over  a  given  velocity 

range  (-1500  to  1500  km/s). 

Several  types  of  corrections  to  the  data  are  necessary  before  actual  line 
strength  measurements  and  profile  modeling  can  be  undertaken.  Corrections 
include:  reseaux  corrections,  cosmic  ray  hit  corrections,  coaddition  of 

multiple  spectra  of  a  common  star,  correction  to  the  flux  for  scattered  light 
between  adjacent  orders,  and  the  fitting  of  a  background  and  stellar 
continuum.  Reseaux  arc  markings  placed  in  the  optics  of  the  telescope  in 
order  to  enable  the  geometric  correction  of  the  spectra,  by  determining  the 
distortion  of  the  optics  due  to  temperature  changes.  The  geometric  correction 
is  done  automatically  at  the  lUE  ground  station  as  the  spectra  arc  reduced  from 
the  satellite  readout.  The  reseaux  markings,  however,  block  sections  of  the 
spectra  (Figure  16)  and  arc  also  easily  seen  as  the  evenly  spaced  white  spots 
on  the  photowritc.  This  "bad"  data  is  flagged  with  a  dotted  line  so  that  it  is  not 
included  in  the  analysis  of  the  absorption  line. 


Figure  16;  A  rescaux  mark  (dotted  line  segment)  in  the  optics  obscured  part  of  the 
spectrum  of  C  11  at  1334  Angstroms. 

Occasionally  the  "bad"  data  quality  flag  extends  to  segments  of  the 
spectra  which  arc  in  fact  good  data  as  can  be  determined  from  comparison 
with  other  spectra  toward  the  same  star  or  through  simple  visual  inspection. 
(Figure  17).  The  improper  data  quality  flags  must  be  changed  in  the  spectral 
line  nic. 
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Figure  17:  An  improperly  flagged  region  (shown  as  the  dotted  segment)  of  the  A1  III  line 
at  1862  Angstroms  which  is  actually  good  data. 
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Conversely,  regions  of  a  line  where  a  cosmic  ray  has  hit  the  detector,  resulting 
in  a  noise  spike,  have  to  be  reflagged  as  bad  data  (Figure  18)  and  these  ray  hits 
appear  as  very  black  spots  on  the  photowrite  image. 


Figure  18;  The  enormous  flux  spike  found  at  the  location  of  a  cosmic  ray  hit  on  the 
detector. 

Coaddition  (the  averaging  of  the  flux  values  for  the  two  lines)  of 
multiple  images  of  the  same  star  has  several  advantages.  The  two  observations 
of  HD  123884,  for  instance,  were  taken  with  an  intentional  velocity  offset  so 
that  the  rescaux  marks  would  not  obscure  the  same  part  of  the  spectrum  in 
both  observations.  Through  coaddition,  the  signal  to  noise  ratio  improves  by  a 
factor  of  fl  and  most  of  the  spectrum  will  have  "good"  data  quality  flags.  In 

order  to  coadd  spectra,  the  velocity  offset  must  be  corrected  so  that  both 
spectral  lines  arc  centered  at  the  proper  velocity. 


The  lUE  spectrum  consists  of  many  orders  which  overlap  slightly  in 
wavelength  at  the  ends  of  adjacent  orders.  The  proximity  of  the  orders, 
especially  in  the  shorter  wavelength  region,  causes  the  scattering  of  light 
from  one  order  onto  the  order  next  to  it.  The  scattered  light  gives  artificially 
high  flux  for  which  the  software  compensates  by  setting  the  zero  flux  level  too 
high.  The  scattered  light  correction  is  determined  by  looking  at  known 
saturated  lines  which  should  have  zero  flux  at  their  bottom  and  shifting  the 
zero  flux  point  to  this  level.  We  apply  the  correction  to  all  of  the  lines  in  that 
order  and  interpolate  for  orders  containing  no  saturated  lines.  The  scattered 
light  correction  introduces  uncertainty  in  the  data  because  saturated  lines 
have  noise  and  arc  not  perfectly  flat  at  the  zero  flux,  so  the  level  must  be 
estimated  through  visual  inspection  or  through  taking  an  average  of  the  flux 
over  the  "flat"  region.  Many  orders  have  no  saturated  lines  yet  arc  obviously 
affected  by  scattered  light  and  must  be  corrected  with  the  best  approximation 
(Figure  19). 


Figure  19  :  An  example  of  a  line  in  an  order  with  scattered  light.  The  1260  Angstrom  Si  II 
line  is  saturated,  and  the  bottom  of  the  lines  serves  as  a  reference  for  the  flux  correction. 
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Before  a  spectrum  can  be  analyzed  for  interstellar  features,  the 
background  continuum  of  light  from  the  star  must  be  fit  to  the  spectral  line. 
The  spectrum  is  then  normalized  by  the  background  so  that  the  light  from  the 
star  has  a  flux  value  of  1  and  regions  where  light  has  been  absorbed  by 
interstellar  gas  have  values  less  than  1.  with  0  being  the  maximum  absorption 
for  saturated  lines  in  which  all  of  the  light  at  a  given  wavelength  has  been 
absorbed  (Figure  20). 


Figure  20;  The  continuum  fit  for  the  Si  II  1526  line  in  HDl  19608  and  the  normalized 
srecirum.  The  boxes  show  the  locations  of  cubic  spline  and  the  dotted  line  indicates  the 
region  excluded  fron.  the  fit. 
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The  procedure  lo  set  the  continuum  fit  takes  a  given  number  of  cubic 
splines  and  fits  a  continuum  based  on  entered  parameters  which  are:  velocity 
range  of  fit,  number  of  points  over  which  to  average  for  each  spline,  number 
of  cubic  splines  to  be  calculated  for  the  fit,  regions  of  the  line  to  be  excluded 
from  calculations,  the  number  of  iterations  of  the  fit,  and  whether  the  splines 
are  calculated  at  evenly  spaced  intervals  or  include  the  edges  of  excluded 
regions.  For  any  given  spectrum,  there  can  be  many  possible  fits  to  the 
continuum.  In  regions  where  the  background  continuum  is  relatively  even,  a 
fit  is  simple  and  the  uncertainty  introduced  by  the  fit  may  be  only  a  few 
percent.  On  the  other  hand,  in  regions  where  the  stellar  profile  shows 
absorption  from  the  atmosphere  of  the  star,  a  unique  fit  is  difficult  or 
impossible,  introducing  considerable  uncertainty  in  the  continuum  fit  and  the 
normalized  spectral  line.  Results  of  the  uncertainty  include  systematic  errors 
in  column  density  caused  by  an  offset  or  incorrect  central  velocity  caused  by 
an  incorrect  slope.  The  distinction  between  the  stellar  and  the  interstellar 
absorption  is  a  complicated  question  worthy  of  a  separate  section  of  discussion. 

Stellar  Ultraviolet  Analysis 

The  background  stars  observed  for  this  investigation  in  the  ultraviolet 
were  hot  distant  stars  ranging  in  spectral  type  from  B1  lo  B8/9.  With  surface 
temperatures  from  10,000  to  25,000  Kelvin,  these  stars  emit  strong  ultraviolet 
radiation  making  them  excellent  background  sources  for  the  study  of 
interstellar  gas,  and  are  luminous  enough  to  be  observed  at  long  distances. 

One  difficulty  present  in  the  study  of  interstellar  absorption  lines,  however,  is 
ih.'.t  the  atmospheres  of  the  stars  thcm.selvcs  also  ab.sorb  certain  wavelengths 
in  the  ultraviolet.  Accurate  analysis  of  the  abundances  of  ions  in  interstellar 
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lines  requires  that  no  contamination  from  the  star's  own  ions  be  included  in 
the  interstellar  calculations.  At  this  time,  work  is  in  progress  to  establish 
standard  ultraviolet  classifications  of  spectral  type  (Rountree  and  Sonneborn 
1991)  similar  to  the  standards  which  have  been  established  in  the  optical 
portion  of  the  spectrum  (Morgan  et  al.  1943).  Some  of  the  spectral  types  have 
been  well  observed  and  have  been  designated  ultraviolet  standard  stars, 
especially  those  B  type  stars  of  average  luminosity.  Other  less  well  observed 
superluminous  type  stars  have  yet  to  be  set  to  definite  standards.  When 
standards  are  available,  one  can  compare  the  standard  star  of  a  given  spectral 
type  to  an  observed  star  of  the  same  type  and  the  stellar  lines  should  be  the 
same.  Since  the  standard  spectra  should  be  taken  of  nearby  stars  with 
virtually  no  interstellar  lines,  the  interstellar  absorption  lines  of  an 

observed  star  would  be  obvious  in  comparison  with  the  standard. 

The  problem  is,  in  fact,  more  complex  as  the  shape  of  stellar  absorption 
lines  depends  on  many  different  factors  of  the  star's  dynamics.  The  rotation  of 

the  star  broadens  the  lines  and  generally  makes  them  symmetrical.  The 
velocity  of  the  star  toward  or  away  from  us  causes  a  Doppler  shift  of  all  of  the 
stellar  lines  by  the  same  amount  corresponding  to  the  star's  radial  velocity. 

The  amount  of  broadening  and  the  shape  of  stellar  lines  depends  on  where  in 
the  atmosphere  the  ions  are  found;  the  ions  at  the  outer  edges  of  the  star's 
atmosphere  tend  to  give  th»''  broadest  lines  as  this  gas  rotates  the  fastest. 
However,  spectra  of  species  in  the  outer  atmospheres  of  some  stars  are  subject 
to  strange  shapes  caused  by  the  asymmetric  outflow  of  particles  from  some 
stars  which  is  known  as  stellar  wind. 

Fortunately,  the  strong  stellar  lines  in  the  stars  we  observed  were  very 
broad  and  sloping  compared  to  the  narrower  interstellar  lines.  Strong  stellar 


Silicon  III  lines  were  observed  in  the  spectra  of  all  three  stars.  If  there  is 
interstellar  Si  III,  its  absorption  spectrum  is  completely  dominated  by  the 
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stellar  line.  Excited  fine  structure  Si  II  lines  were  also  seen  which  have 
the  same  shape  as  the  stellar  Si  III  and  look  very  different  from  the  ground 
state  Si  II  interstellar  absorption  (Figure  21).  We  know  that  the  excited  fine 
structure  must  be  stellar  because  the  density  and  pressure  of  gas  in 
interstellar  space  would  not  be  enough  to  produce  these  lines  at  a  detectable 
strength;  they  could  only  be  formed  in  a  star. 


stellar  lines  in  HD  123884 
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Figure  21  :  Stellar  Si  III  and  Si  II  absorption  lines  (which  ha^^e  not  been  normalized  since 
the  stellar  light  is  by  definition  the  background). 
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In  many  cases,  an  obvious,  narrow  interstellar  absorption  line 
protrudes  from  the  bottom  of  a  broad,  sloping  stellar  absorption  line.  The 
Carbon  IF  lines  at  1334  and  1335  A  ir  HDl  19608  and  the  Aluminum  Ill  lines  at 
1854  and  1862  A  in  HD123884  are  the  most  notable  examples  (Figure  22). 


Figure  22:  The  unnormalized  absorption  lines  of  C  11  in  HDl  19608  and  A1  III  in  HD123884 
showing  distinguishable  stellar  and  interstellar  components. 

The  most  difficult  lines  to  distinguish  are  those  in  which  stellar  and 
interstellar  components  both  contribute  significantly  to  the  line.  There  is 
almost  certainly  interstellar  A1  III  toward  HDl  19608  and  BD+2  2711  as  well  as  Si 
IV  toward  BD+2  2711,  yet  we  could  only  set  limits  on  the  interstellar 
abundances  because  we  could  not  distinguish  the  stellar  and  interstellar 
components.  Also,  weak  lines  without  a  clearly  stellar  or  clearly  interstellar 
shape  have  generally  not  been  as  well  studied  as  the  strong  lines  that 
disiinguisli  one  spectral  type  from  another.  But,  without  particular 
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characteristics  of  the  physical  properties  associated  with  the  line,  we  cannot 
absolutely  call  such  lines  interstellar.  For  example,  weak  absorption  lines 
with  a  very  small  but  symmetric  and  sloping  shape  are  seen  at  the  positions  of 
2063,  2027,  and  2022  A  (which  arc  the  locations  of  excited  fine  structure  Cobalt 
II  lines)  that  should  not  be  seen  in  the  pressure  and  density  regimes  of  halo 
gas  and  could  therefore  be  small,  weak,  stellar  lines.  A  more  flat  bottomed, 
ground  state  Co  II  line  appears  at  2058  A  which  could  be  stellar,  interstellar,  or 
both;  it  is  not  possible  to  determine  at  this  time  with  the  low  velocity 
resolution  of  the  data.  Higher  resolution  spectra  from  the  Hubble  Space 
Telescope's  ultraviolet  spectrograph  may  help  answer  some  of  these  questions 
in  the  future.  Also,  one  project  that  might  be  extremely  useful  for  both 
astronomers  studying  the  stars  and  those  studying  the  interstellar  gas  would 
be  a  large  scale  study  of  the  characteristics  of  the  respective  line  types  to 
establish  a  standard  model  of  the  stellar  and  interstellar  lines.  Their 
respective  contributions  to  a  given  spectrum  could  then  be  systematically 
identified  and  separated. 

Analysis  of  the  Interstellar  Absorption  Lines 

We  made  the  necessary  corrections  and  conducted  preliminary  analysis 
for  a  number  of  the  clear  interstellar  absorption  lines  of  neutral,  singly 
ionized  and  highly  ionized  gas.  Figure  23  shows  the  strongest  and  most  easily 
detected  species  in  each  star.  Because  of  the  faintness  of  the  distant  star 
HD123884,  two  full  night  exposures,  which  were  then  coadded,  are  still 
underexposed  and  too  noisy  for  accurate  analysis  in  the  shorter  wavelengths 
of  the  ultraviolet  spectrum.  The  long  wavelength  portion  of  the  spectrum. 


Figure  23:  The  strong  ultraviolet  absorption  lines  toward  HDl  19608,  HD123884 

and  BD+2  2711. 
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Figure  23  continued. 


however,  is  well  exposed  and  has  provided  important  information  about  this 
line  of  sight.  The  original  spectra  of  HDl  19608  from  the  archives  had  been 
extracted  from  the  raw  satellite  data  at  Goddard  Space  Center  when  older 
software  systems  were  in  use.  At  our  request,  NASA  re-extracted  the  spectra  of 
HDl  19608  with  the  modem  technique  to  make  its  resolution  comparable  to  our 
other  observed  data.  Upon  receipt  of  the  re-extracted  spectra  in  April  1992,  we 
reanalyzed  most  of  the  lines  toward  HDl  19608  at  the  improved  resolution. 


Ultraviolet  Interstellar  Analysis  with  the  Curve  of  Growth 


Once  the  spectra  have  been  corrected  and  the  continua  have  been  fit 
and  normalized,  a  curve  of  growth  fitting  technique  is  used  to  determine  the 
abundances  of  the  various  ions.  The  first  step  in  the  process  is  to  measure  the 
area  of  the  absorption  line.  Since  the  vertical  height  of  the  line  is  normalized 
to  a  background  of  one,  the  area  measured  has  units  of  velocity  which  arc 


then  converted  to  wavelength  by  scaling  the  value  by  the  factor  (X  /  C)  where 
X  is  the  wavelength  of  the  line  in  Angstroms  and  C  is  the  speed  of  light  in 
km/s.  The  resulting  value  is  called  the  equivalent  width.  The  information 
which  we  want  to  know,  however,  is  not  the  equivalent  width,  which  is  a 
measure  of  line  strength,  but  the  actual  abundance  of  the  gas  in  number  of 
atoms,  or  column  density.  The  relationship  between  the  equivalent  width  and 

the  column  density  is  called  the  curve  of  growth.  In  general,  the  observed 
strength  of  a  line  increases  as  the  column  density  "grows".  For  weak  lines,  the 
strength  of  the  line  is  directly  proportional  to  the  column  density.  As 
absorption  lines  become  saturated  (all  of  the  light  is  absorbed  in  a  region),  the 
curve  becomes  nearly  flat  and  a  large  change  in  the  number  of  atoms  will 
result  in  very  little  change  in  the  line  strength.  Eventually,  if  a  line  is 
extremely  saturated,  then  the  natural  Doppler  broadening  caused  by  the 
average  kinetic  motions  of  atoms  in  the  gas  will  cause  the  formation  of 
"wings"  or  extended  regions  of  the  absorption  line.  The  damping  wings 
increase  in  area  by  a  factor  of  the  square  root  of  the  column  density.  We  used 
the  Macintosh  Excel  and  Crickclgraph  spreadsheet  programs  to  plot  theoretical 
curves  of  growth  for  different  cloud  temperatures,  or  Doppler  broadening  (b). 
Appendix  C  contains  a  more  detailed  mathematical  description  of  the 
theoretical  curves  of  growth.  The  theoretical  curves  are  plots  of  the  log  of  the 
equivalent  width  divided  by  the  wavelength  versus  the  log  of  the  column 
density  multiplied  by  the  wavelength  and  the  oscillator  strength  (a  measure  of 
the  probability  of  an  atomic  transition).  However,  observed  data  is  plotted  as 
ilie  log  of  equivalent  width  divided  by  the  wavelength  versus  the  log  of  the 
wavelength  times  the  oscillator  strength.  The  horizontal  shift  required  to  fit 
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ihc  observed  sel  of  points  for  each  ion  to  the  theoretical  curve  of  growth  is 
then  the  log  of  the  column  density  of  that  ion  (Figure  24). 


log  (  Nj  fjic  X  ) 


Figure  24:  Theoretical  curves  of  growth  at  different  cloud  temperatures,  or  Doppler 
broadening  (b).  The  dark  plotted  points  are  the  observed  equivalent  widths  of  the  line  of  S 
II  toward  BD+2  2711.  The  shift  of  the  data  to  the  theoretical  curve  equals  log  (column 
density  (N)  ). 


This  curve  of  growth  technique  gives  the  total  column  density  of  each  fitted 
ion  along  the  line  of  sight  toward  each  of  the  stars  observed  in  the  ultraviolet. 
Appendix  C  also  includes  all  of  the  fitted  curves  of  growth. 

Table  2:  Column  Densities  by  Curve  of  Growth  Fit 


Ion 

BD  +2  2711 

HD119608 

HD123884 

Si  II 

3.55  E  14 

1.00  E  15 

2.00  E  15 

S  II 

1.00  E  15 

5.01  E  15 

Zn  II 

3.55  E  12 

1.78  E  13 

2.00  E  13 

CrII 

1.26  E  13 

5.01  E  13 

3.98  E  13 

Ni  II 

6.31  E  13 

CIV 

6.31  E  13 

8.91  E  13 

Si  IV 

2.29  E  13 

N  V 

5.37  E  13 

The  total  column  density  is  the  number  of  absorbing  atoms  along  the  line  of  sight  in  units 
of  number  of  aioms/cm^ .  The  main  source  of  uncertainty  of  the  column  density  is  the 
uncertainty  of  equivalent  width  typically  20  to  30%. 

The  principal  experimental  uncertainties  in  the  column  densities  lie  in 
the  uncertainty  of  the  equivalent  width  caused  primarily  by  the  uncertainty 
of  the  continuum  fit  which  can  cause  about  20%  variation  in  the  measured 
equivalent  width  of  strong  lines.  The  uncertainty  of  column  density  is 
proportional  to  the  uncertainty  in  equivalent  width.  in  addition,  there  is  some 
uncertainty  in  the  best  fit  of  the  curve  of  growth,  particularly  with  species  in 
the  "flat",  or  saturated  portion  of  the  curve.  To  help  minimize  uncertainty, 
with  the  exception  of  the  single  N  V  line,  we  only  fit  ions  that  have  multiple 
absorption  lines.  The  average  uncertainty  in  fitting  a  set  of  points  is  much 
less  than  that  of  a  single  data  value. 
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Spectral  Analysis  with  Profile  Fitting 

The  curve  of  growth  fitting  technique  gives  the  overall  column  density 
of  each  ion,  along  the  line  of  sight,  but  to  determine  the  column  density  of 
each  cloud  component  by  velocity,  we  again  used  a  profile  fitting  technique. 
Toward  HD  119608,  the  cloud  model  derived  from  the  high  resolution  optical 
data  (5  km/s  velocity  resolution)  can  be  applied  to  the  low  resolution  (  25 
km/s)  ultraviolet  data.  We  varied  the  column  density  and  added  additional 
cloud  components  at  higher  velocities  in  order  to  fit  the  profiles  of  several  of 
the  clear  and  interesting  interstellar  lines.  An  example  is  given  in  Figure  25. 

The  profile  fitting  process  gives  the  individual  velocity,  Doppler 
broadening  width,  and  column  density  of  each  separate  cloud  along  the  line  of 
sight.  The  resulting  cloud  column  densities  are  listed  in  Table  3.  The 
uncertainties  in  these  values  are  proportional  to  the  uncertainties  in  the 
equivalent  widths,  typically  20  to  30%  for  these  data. 


SWP5536HLC+SWP5647HSM 


Cloud  Model 


Figure  25;  Profile  fit  of  the  Si  II  lines  toward  HDn9608  using  the  cloud  model  from  the 
optical  as  the  starting  model  for  the  fitting  process.  The  model  profiles  are  shown  with 
the  observed  lines. 


48 


Table  3: 


Column  Densities  by  Profile  Fitting 


C4  K>K _ C»  H-W _ _ T»t» 


SUf  CpmpwMSt 

(kma-*)  (iO 

H 

'  «n“ 

b 

*J  (km** 

•ifii 
*)  (km*-* 

S 

)  (IO”  cm’* 

b 

(km  •' 

)  (kmi")  (10 

N 

'  cm’* 

b 

(km  •’* 

bmt 

(I0“  cm-*) 

KO  I19608  1 

-39.3 

03 

3.0 

-39  3 

03 

30 

-33.3 

07 

40 

-33  3 

04 

40 

-33  3 

10 

40 

-188 

1.4 

40 

-16  8 

3  1 

4  0 

-16  0 

14 

40 

-7.3 

40 

40 

-7.4 

44 

40 

-74 

30 

40 

3.3 

88 

30 

30 

7.4 

SO 

37 

87 

44 

10  7 

0.8 

30 

tOT 

11 

3  0 

107 

03 

30 

19.7 

94 

80 

19  7 

93 

4.3 

31  0 

10  4 

7  1 

31.7 

30 

4.0 

304 

31 

34 

34.4 

14 

40 

IID  UTt80  3 

-16  8 

04 

40 

-S7 

4  0 

40 

33 

14 

30 

10  7 

10 

40 

110  118346  3 

-18  7 

3  1 

38 

30 

-4  7 

33 

4  4 

— 

39  3 

1  0 

40 

HO  130086  1 

-38  9 

07 

60 

-37  3 

09 

60 

14 

-184 

4  7 

30 

-17  7 

44 

4.4 

-14  3 

47 

80 

-4  9 

7 

4  7 

-4  7 

7  1 

38 

-17 

39 

SO 

3  1 

08 

60 

33 

1  3 

40 

14  1 

1 

80 

<0  7 

8.0 

KO  133884  1 

-39.3 

0  4 

30 

<06 

30 

-33  3 

09 

40 

<1  0 

4  0 

-188 

30 

40 

-16  8 

0  9 

4  0 

-4  7 

9.0 

40 

-6  8 

too 

40 

33 

30 

30 

3.3 

30 

30 

10.7 

40 

4.0 

10  7 

40 

4.0 

31.3 

140 

7.0 

31.4 

33.0 

7.0 

31.7 

4.0 

8.0 

31.7 

30 

80 

Uir 

Ion 

Cay*oan«(u 

vlSft 

0 

N 

mm 

m 

HOl 19600 

Si  II 
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•600 

IIHI 

mm 

2 

24 

20813 

3 

•360 

6 

6Qri3 

'29  i 

2 

60ri3 

5 

•?5  3 

4 

aori3 

•<66 

4 

I.6ri4 

*7  3 

S 

«98I« 

22 

3 

60rl4 

107 

2 

2.0FI4 

19  7 

6 

22714 

31  7 

80712 

26 

1  3713 

HOl 19600 

, 

'640 

3 

2 

•140 

99 

90ri  1 

5 

•29  3 

2 

70fi2 

4 

•25  3 

u 

•1 

9  VI2 

'  •• 

i 

’  S?I2 

27 

5 

70712 

2 

■i07l2 

c 

30f»2 

.1 

2  >«2 

•’ 

J 

•  0r»2 

HOl 19600 

Cf  II 

1 

-J9C 

j 

2  0712 

2 

•56C 

V 

2  371? 

i 

-29  5 

2 

2a»l7 

2271? 

d 

•IC^I2 

6 

3 

6?^l2 

70712 

»0  7 

7 

6  7712 

6 

99712 

MB 

.1 

•:CV>2 

In  order  to  compare  the  precision  of  the  two  methods  for  obtaining 
column  density,  Table  4  gives  the  sum  of  the  column  densities  of  the  cloud 
models  and  the  total  column  density  derived  from  the  curve  of  growth 
iccliniquc.  Tliesc  values  are  comparable  to  within  their  uncertainties. 
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Table  4:  Total  Column  Densities  by  Two  Methods 


(N  in  atoms  /  cm^) 


HDl 19606 

Ion 

(Prof  lie  Fit) 

tv' 

(Curve  of  Growth) 

N 

^Diff  erence 

Average 

Si  II 

I.87el5 

1.00  el  5 

47 

I.43ei5 

Zn  II 

3.61  613 

1.78  el  3 

50 

2.73  el  3 

■  Cr  II 

4.25  el  3 

5.01  el  3 

1  6 

4.63  el  3 

VI  Results 

Taken  together,  the  radio,  optical  and  ultraviolet  observations  of  this 
unusual  region  in  our  Milky  Way  galaxy,  lead  to  four  principal  results,  all  of 
which  are  suggestive  of  a  large  scale  flow  dynamic  such  as  a  galactic  chimney. 
First,  we  see  a  large  scale  outflowing  cloud  in  neutral  hydrogen  emission 
spectra  spanning  a  distinct  region  of  the  sky  across  the  area  from  310®  to  330® 
galactic  longitude  and  from  about  38®  to  48®  galactic  latitude.  After  we 
identified  the  only  region  of  marked  outflow  in  the  northern  galactic  sky  in 
the  low  resolution  Bell  Labs  H  I  survey,  wc  mapped  the  location  at  high  spalia. 
resolution  using  the  140  foot  National  Radio  Astronomy  Observatory  radio 
telescope.  The  high  resolution  mapping  confirms  the  large  scale  outflow  of 
neutral  hydrogen  and  shows  fine  scale  structure  in  the  motions  of  interstellar 
gas.  We  determine  the  distance  limits  of  the  outflow  using  high  resolution 
optical  absorption  spectra  of  stars  along  the  line  of  sight  through  the 
outflowing  hydrogen  region.  The  positive  velocity  components  in  the  line 
profile  models  of  singly  ioni/cd  calcium  appear  at  a  height  of  about  1  kpc 


50 

above  the  galactic  plane  and  increase  in  strength  into  the  halo.  The  column 


density  of  the  outflowing  components  increases  from  the  height  3.1  kpc  at 
HDl  19608  to  the  star  HD123884,  believed  to  be  at  8.9  kpc  above  the  plane  of  the 
galaxy.  The  relation  between  height  and  outflowing  column  density  can  be 
seen  in  Table  5.  Although  it  is  unlikely  that  the  cloud  extends 


Table  5:  Outflow  Column  Density  with  Distance 


Star 

Distance 

Ca  II  N  at  23km/s 

1 17880 

0.33  kpc 

0 

1 18246 

1.1  kpc 

l.Oel  1 

1 19608 

3.1  kpc 

1 1.3  el  1 

123884 

8.9  kpc 

20.0  e  1 1 

that  far  above  HDl  19608,  this  evidence  certainly  indicates  that  gas  is  flowing 
out  of  the  plane  of  the  galaxy  to  heights  thousands  of  light  years  beyond  the 
disk.  The  detections  of  the  outflow  feature  in  singly  ionized  titanium  confirms 
that  the  observations  of  Ca  11  and  H  I  indeed  represent  the  same  neutral  cloud. 
Interestingly,  we  do  not  see  the  strong  outflow  of  H  I  toward  the  star  HD123884 
which  is  observed  in  the  calcium  absorption  toward  this  star.  HD123884  may 
mark  the  edge  of  a  line  of  .sight  where  the  accelerated  cool  cloud  becomes 
ionized. 

In  addition  to  the  large  scale  regional  outflow  of  the  gas  elouds  in  the 
HDl  19608  line  of  sight,  the  .second  principal  result  of  this  study  i.s  the  detection 


of  certain  ions  in  the  ultraviolet  spectra  which  provide  some  crucial  keys  to 
unravelling  the  mystery  of  this  galactic  flow.  Highly  ionized  species  detected 
in  ultraviolet  absorption  generally  indicate  hot  gas.  We  see  strong  C  IV  in  the 
spectra  toward  the  stars  HDl  19608  and  BD+2  2711  and  also  delect  the  high  ion  Si 
IV.  The  approximate  C  IV  /  Si  IV  ratio  from  the  curve  of  growth  values  of 
column  density  is  3.9  toward  HDl  19608.  While  the  C  IV  and  Si  IV  are  not 
unique  to  this  line  of  sight,  the  most  highly  ionized  species  in  the  absorption 
spectrum  of  HDl  19608  is  the  ion  N  V  which  has  not  been  found  in  the 
interstellar  absorption  spectra  of  any  other  high  latitude  star.  Nitrogen  V  may 
be  a  particularly  important  indicator  of  gas  temperature  as  it  is  generally 
believed  to  indicate  an  area  of  collisionally  ionized  gas  cooling  through  a 
narrow  temperature  range  of  about  10^  Kelvin. 

While  the  presence  of  high  ions  is  a  significant  result  in  and  of  itself, 
just  as  interesting  is  the  third  result:  velocity  trends  in  motion  that  we  observe 
in  the  high  ions.  HDl  19608,  HD125924  and  HD121968  all  show  a  trend  of 
stronger  absorption  at  infalling  velocities  for  the  high  ions  and  outflowing 
for  the  low  ions.  This  trend  is  not  observed  in  the  higher  latitude  star  BD+2 
2711,  where  the  absorption  lines  are  centered  at  zero  velocity.  Therefore,  the 
flow  of  cold  gas  out  of  the  disk  of  the  galaxy  and  hotter  gas  into  the  disk  occurs 
up  to  the  latitude  of  the  former  stars  (43  to  56')  but  docs  not  extend  to  the 
height  of  BD+2  2711  at  62'. 

These  flows  are  large-scale  phenomena,  but  they  arc  localized  in  a  well 
defined  region  of  space.  The  observations  are  consistent  with  a  chimney 
scenario  and  provide  important  constraints  for  future  theoretical  models.  Gas 
flowing  out  of  the  disk  in  a  galactic  chimney  is  expected  to  be  hot,  but  models 
(eg  Norman  and  Ikeuchi  1989)  cannot  confidently  predict  an  exact 
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temperature  range.  Our  observation  of  the  high  ions  C  IV,  S  IV,  and  especially 
N  V  flowing  toward  the  disk  suggests  that  the  gas  was  hotter  yet,  perhaps  10^ 
Kelvin,  when  originally  ejected  upwards.  The  outflowing  neutral  cloud  that 
we  observe  is  therefore  not  cold  condensations  within  that  chimney  flow  but 
could  be  a  supernova  shell  or  swept  up  ambient  disk  gas  propelled  upwards  by 
a  later  chimney  event. 

The  fourth  principal  result  of  this  study  is  the  strong  abundance  of 
singly  ionized  chromium  in  addition  to  unusually  strong  calcium  and  titanium 
absorption  along  this  line  of  sight.  All  of  these  refractory  elements  condense 
very  easily  onto  interstellar  grains  of  dust  leaving  few  atoms  of  the  element  in 
the  gas,  so  the  absorption  line  is  usually  very  weak  and  indicates  only  a 
fraction  of  the  actual  number  of  atoms  in  space  because  most  are  condensed 
onto  grains.  When  the  grains  have  been  destroyed  by  some  mechanism  such 
as  the  shock  from  a  supernova  or  high  velocity  clouds  (Routly  and  Spitzer 
1952)  associated  with  a  chimney  or  fountain  flow,  the  normally  "depleted" 
elements  are  liberated  into  the  gas,  or  "dedepleted",  and  can  therefore  cause  a 
strong  absorption  line.  The  pattern  of  chromium  depletion  is  similar  to  that  of 
the  optically  observed  Ca  II  and  Ti  II.  Like  Ti  II,  Cr  II  is  a  dominant  ionization 
stage  in  regions  of  neutral  hydrogen.  From  the  column  densities  in  the 
profile  fitting  technique  we  find  that  the  ratio  of  the  abundance  of  chromium 
varies  from  3  to  84%  of  the  solar  abundance.  Chromium  is  panicularly 
enhanced  in  both  the  outflowing  and  infalling  velocity  components,  again 
consistent  with  grain  disruptions  in  infalling  fountain  material  or  in 
supernova  shocks  (Table  6). 
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Table  6  :  Abundances  of  Refractory  Elements 

comnonent  eas  phase  abundance.  N(x)/N(H  1) 

(LSR  velocity)  relative  to  the  solar  abundances 


Cr  II  Ca  II  Ti  II 


-29.3 

0.84 

0.002 

-233 

0.52 

0.003 

0.0130 

-16.8 

0.17 

0.002 

0.056 

-7.3 

0.03 

0.0005 

0.018 

2.2 

0.07 

0.001 

0.046 

10.7 

0.34 

0.0008 

0.017 

19.7 

0.15 

0.003 

0.J40 

31.7 

0.49 

0.002 

0.106 

Unlike  chromium,  zinc  does  not  condense  easily  onto  grains  and 
therefore  serves  as  a  good  indicator  of  actual  metal  content  of  the  gas  (Pettini 
ct  al.  1990).  Since  the  ionization  potential  of  Zn  II  to  Zn  III  is  17.96  eV,  which  is 
much  greater  than  the  13.6  eV  needed  to  ionize  hydrogen,  Zn  II  is  also  the 
dominant  ionization  stage  in  neutral  hydrogen  regions.  The  measured  ratio  of 
zinc  with  respect  to  hydrogen,  Zn  II  /  H  I,  is  essentially  its  total  abundance, 
Zn/H.  Surveys  of  interstellar  zinc  to  hydrogen  ratios  show  the  average  value 
to  be  2.61x10'^,  or  58%  of  the  solar  abundance  (Harris  ct  al.  1983).  The  total 
column  density  of  zinc  in  the  region  of  this  study  is  average  but  the  individual 
column  density  at  the  positive  velocity  clouds  is  158%  of  solar  and  the  negative 
velocity  component  is  84%  of  solar,  while  the  zero  velocity  material  has  a  zinc 
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abundance  of  only  23%  of  solar.  This  work  is  the  first  study  of  interstellar  zinc 
abundances  in  individual  cloud  components  and  reveals  that  the  zinc 
abundance  is  more  complex  than  observed  in  previous  general  surveys.  At 
least  along  this  line  of  sight,  the  abundance  varies  significantly  with  velocity 
and  may  correlate  with  different  histories  of  various  individual  clouds.  Both 
the  outflowing  and  infalling  clouds  in  this  region  show  particularly  enhanced 
zinc,  and  thus  high  metallicity,  characteristic  of  gas  recently  cycled  through  a 
supernova. 

Of  course,  these  results  depend  upon  the  continuum  fitting  process,  the 
resolution  of  the  absorption  lines  and  the  reliability  of  the  profile  fitting 

method.  However,  even  if  the  zinc  column  density  were  reduced  by  a  factor  of 
two,  the  basic  conclusions  would  still  stand.  In  this  case,  the  average  zinc 
abundance  in  the  positive  velocity  components  would  be  84%  of  the  solar 

abundance,  which  is  still  a  factor  of  1.4  greater  than  the  average  interstellar 
value. 

With  the  results  of  a  strong  outflowing  feature  with  spatial  and  distance 
constraints,  presence  of  high  ions  including  N  V,  a  velocity  separation 

between  high  and  low  ions,  and  enhanced  chromium  and  zinc  at  the  positive 

velocity  components,  we  have  accomplished  our  preliminary  objectives  of 
exploring  the  candidate  chimney  region  in  space,  distance,  velocity  and 


ionization. 


VII  Conclusions 


No  static  halo  model  can  explain  our  very  dynamic  observations.  Strong 
activity  is  accelerating  material  causing  a  separation  in  the  velocity  of  the  low 
and  high  ions.  Definite  trends  have  been  noted  in  space,  distance,  velocity, 
and  ionization.  Enhanced  chromium  and  calcium  indicated  grain  destruction 
associated  with  a  shocked  region  while  high  zinc  abundance  reveals  high 
metallicity  characteristic  of  gas  recently  recycled  through  a  supernova.  The 
strong  feature  of  outflowing  neutral  hydrogen  and  low  ions  is  consistent  with 
the  swept  up  shell  of  a  supernova  remnant  or  could  be  an  unrelated  feature. 
The  ultraviolet  results  are  consistent  with  a  very  hot  chimney  model  with  a 
very  long  cooling  time  in  which  the  gas  would,  in  fact,  have  to  be  too  hot  to 
observe  with  our  methods  on  the  way  out  of  the  disk  (perhaps  one  million 
degrees  Kelvin)  and  we  sec  the  returning  fountain  flow  only  when  the  gas  has 
cooled  to  the  relatively  cooler  "high"  ions  that  we  see  at  predominantly 
infalling  velocities.  If  this  local  phenomenon  is  typical  of  the  galaxy  as  a 
whole,  then  the  cycling  of  gas  between  the  stars  in  the  disk  and  the  halo  must 
play  an  important  role  in  the  dynamics  and  evolution  of  the  Milky  Way. 
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Appendix  A 

Neutral  Hydrogen  Emission  Contour  Maps 
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Appendix  B 

Atomic  Line  Data  File 


Line 

Wavelength 

Oscillator  Slrenth  (F) 

Upper 

Lower 

Order 

6 

4 

115 

M  I 

1199.5496 

0.133000 

3 

1 

115 

SI  III 

1206.500 

1.660000 

3 

1 

114 

SI  III 

1206.500 

1.660000 

4 

2 

112 

N  V 

1238.821 

0.152000 

4 

2 

111 

M  V 

1238.821 

0.152000 

2 

2 

111 

N  V 

1242. 8C4 

0.075700 

2 

4 

110 

S  II 

1250.584 

0.005350 

4 

4 

110 

S  II 

1253.811 

0.010700 

6 

4 

110 

S  II 

1259.519 

0.015900 

6 

4 

109 

S  II 

1259.519 

0.015900 

4 

2 

109 

SI  II 

1260.4221 

0.960000 

0 

0 

107 

Si  III 

1294.540 

0.000000 

0 

0 

106 

Si  III 

1294.540 

0.000000 

3 

5 

106 

0  I 

1302.1685 

0.048600 

2 

2 

106 

Si  II 

1304.3702 

0.147000 

4 

2 

103 

C  II 

1334.5323 

0.118000 

10 

8 

103 

C  II* 

1335.7077 

0.118000 

4 

2 

99 

Si  IV 

1393.755 

0.528000 

2 

2 

98 

Si  IV 

1402.770 

0.262000 

0 

0 

97 

Si  III  - 

1417.243 

0.000000 

2 

2 

90 

Si  II 

1526.7066 

0.230000 

4 

2 

89 

C  IV 

1548.195 

0.194000 

2 

2 

89 

C  IV 

1550.770 

0.097000 

8 

10 

86 

Fe  II 

1608.4511 

0.062000 

3 

1 

83 

Al  II 

1670.7874 

1.880000 

3 

1 

82 

A1  II 

1670.7874 

1.880000 

4 

2 

76 

Si  II 

1808.0126 

0.005500 

4 

2 

74 

Al  III 

1854.7164 

0.539000 

2 

2 

74 

Al  III 

1862.7895 

0.268000 

3 

1 

83 

C  I 

1656.928 

0.136000 

3 

1 

88 

C  I 

1560.310 

0.081000 

3 

1 

108 

C  I 

1280.135 

0.027800 

3 

3 

106 

0  I* 

1304.858 

0.048500 

3 

1 

106 

0  I** 

1306.029 

0.048500 

2 

2 

115 

Si  II 

1193.289 

0.500000 

2 

2 

116 

Si  II 

1193.289 

0.500000 

4 

2 

116 

Si  II 

1190.416 

0.250000 

8 

6 

79 

N  II 

1751.910 

0.040000 

6 

6 

79 

N  II 

1741.549 

0.067900 

6 

6 

81 

N  II 

1709.600 

0.047000 

6 

6 

95 

H  II 

1454.842 

0.051500 

4 

6 

100 

N  II 

1370.132 

0.130000 

4 

6 

101 

N  II 

1370.132 

0.130000 

2 

2 

67 

Zn  II 

2062.005 

0.202000 

3 

1 

108 

C  I 

1277.245 

0.156000 

8 

6 

67 

Cr  II 

2055.596 

0.109000 

4 

2 

68 

Zn  II 

2025.483 

0.412000 

/ 

The  atomic  line  data  file  gives  the  parameters  about  a  given  line  necessary  for  modeling 
The  oscillator  strengths  and  numbers  of  upper  and  lower  transitions  are  from  Morton 
1991.  The  order  is  the  line  on  the  photowrite  that  contains  the  wavelength  of  interest. 
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Appendix  C:  Curve  of  Growth  Theory 

The  relationship  between  the  equivalent  width  and  the  column  density  is  the 
curve  of  growth. 

Equivalent  width  is  defined  as  the  area  under  the  absorption  line  and,  if  the 
background  is  normalized  to  unity,  it  is  expressed  in  units  of  wavelength,  or 
Angstroms.  Mathematically,  equivalent  width  is  written  as: 


where  lv(0)  is  the  background  or  continuum  intensity  and  ly  is  the  intensity 
at  each  point  along  the  absorption  line.  dX.  is  the  small  change  in  wavelength 
across  the  line,  so  the  sum  of  the  normalized  line  depth  (l-Iy/IyfO))  times  the 
increment  in  width  equals  the  area  under  the  line. 


Equivalent  width  relates  to  the  optical  depth,  Ty,  of  the  cloud  by  the  expression: 


where  the  integral  is  now  expressed  as  a  function  of  frequency  rather  than 
wavelength.  The  optical  depth  of  a  line  is  a  measure  of  the  atomic  element's 

ab.sorption  across  a  line  profile  9(v): 


T  —  Nj  s  (p(v)  where  Nj  is  the  number  of  atoms  in  the  lower 

energy  state  along  a  line  of  sight 

s=  2.65x1 0'2  fjk  where  fjk  is  the  oscillator  strength,  a 

measure  of  the  probabity  of  the  transition 


Ir  a  weak  line,  Ty  «  1,  we  can  expand  the  exponent  so  that 


t„dv  =  8,85xl0'^ 


In  the  case  of  the  weak  line,  then,  the  eqiiivaleni  width  over  the  wavelength  is 
directly  proportional  to  the  column  den.siiv  and  the  o.scillator  strength.  This 
is  the  "linear"  portion  of  the  curve  of  growth  and  the  shape  of  the  line  is 
determined  by  Doppler  broadening. 


For  stronger  lines,  the  optical  depth  is  large  and  the  amount  of  energy 
absorbed  is  no  longer  proportional  to  optical  depth.  The  line  is  said  to  be 
saturated  and  the  observed  equivalent  width  does  not  change  significantly 
over  a  wide  range  of  column  density.  This  is  the  "flat"  portion  of  the  curve  of 
growth.  In  this  case,  the  curve  depends  on  the  detailed  line  profile  and  can  be 
written  as  : 


Wj^  =  2  p  F(C) 


where: 


2  2 
V  m 


F(C) 


.j(,. 


dx 


N.S 

C= _ i— 


Finally,  when  the  optical  depth  becomes  very  large,  radiation  damping  wings 
become  important  and  ultimately  dominate  the  line  shape.  This  is  the  "square 
root"  portion  of  the  curve  of  growth  since,  in  this  limit, 

W  =  2(Nj  sy)2 

where  y  ,  a  measure  of  the  width  of  the  upper  level  of  the  transition,  is 
determined  by  the  Einstein  A  coefficients: 


Combined,  these  equations  result  in  the  theoretical  curve-of-growth  shown 
schematically  in  Figure  A.l. 

As  described  in  the  text,  the  shift  between  observed  data,  plotted  as  log  (W;^,  /  X) 
versus  log  (fjk  X)  and  the  theoretical  curve  of  growth  gives  the  column 
densities  of  the  observed  lines.  Table  A.l  lists  the  observed  equivalent  widths 
of  the  ultraviolet  absorption  lines  and  Figure  A. 2  illustrates  the  curve  of 
growth  analysis.  The  resulting  column  densities  were  presented  in  text  (Table 
2). 
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Table  A.l  :  Equivalent  Widths 
Wavelength  W  HD119608  W  BD-t-2  2711  W  HD123R84 


Si  II  1808 

0.140 

0.515 

0.241 

1526 

0.407 

0.353 

0.407 

1304 

0.348 

0.183 

1260 

0.561 

0.406  • 

S  II  1250 

0.185 

0.0059 

1253 

0.279 

0.131 

1259 

0.280 

0.159 

Zn  II  2026 

0.322 

0.0669 

0.260 

2062 

0.183 

0.034 

0.173 

Cr  II  2065 

0.085 

0.027 

0.125 

2062 

0.165 

0.0476 

0.226 

2055 

0,075 

Ni  II  1751 

0.0822 

0.234 

1741 

0.182 

1709 

0.102 

0.351 

1454 

0.039 

CIV  1548 

0.168 

0.375 

15  50 

0.0968 

0.178 

Si  IV  1393 

0.194 

1407 

0.0896 

N  V  1238 

0.1  15 

Figure  A.l:  Schematic  example  of  weak,  saturated,  and  damped  lines. 


Figure  A. 2:  Fitted  Curves  of  Growth 
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